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The QCD thermodynamics on the lattice provides fundamental theoretical grounds to analyze the 
various experimental data in relativistic heavy ion collisions. So far, most of the numerical simu- 
lations on the lattice have been performed by using the staggered-type fermion actions. Therefore 
it is important to carry out studies using different fermion formulations to test the uncertainties of 
the lattice QCD results. For this purpose, we perform systematic simulations of two-flavor QCD 
with an improved Wilson quark action to investigate the equation of state. We report the current 
status of our project and show the preliminary results of the Taylor expansion coefficients of the 
thermodynamic grand partition function in terms of chemical potential. 
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1. Introduction 

Numerical simulation of lattice QCD at non-zero temperature {T) and quark chemical poten- 
tial {fLq) is an essential tool for quantitative understanding of the QCD phase transition. So far, 
most of the studies have been performed using the staggered-type quaik actions with the fourth- 
root trick of the quark determinant. To test the uncertainties of the lattice QCD results due to 
different fermion formulations and to obtain basic information to analyze the experimental data in 
relativistic heavy ion collisions, systematic studies of the QCD thermodynamics using a Wilson- 
type quark action are called for. Such a study at T 7^ and Hq = has been initiated six years ago 
using the Iwasaki (RG) improved gauge action and the Nf = 2 clover improved Wilson quark ac- 
tion by the CP-PACS Collaboration [|l], ^. The phase structure, the transition temperature and the 
equation of state have been investigated in detail, and also the crossover scaUng in the region near 
the chiral phase transition point has been tested. It is confirmed that a subtracted chiral condensate 
satisfies the scaling behavior with the critical exponents and scaling function of the 3-dimensitonal 
0(4) spin model, suggesting the chiral phase transition is in the same universality class as the 0(4) 
spin model. Moreover calculations of various physical quantities at T = such as the light hadron 
masses have been carried out using the same action 

Since there are numbers of technical progresses in treating the system with finite baryon den- 
sity in the past six years, it may be worth while to revisit the QCD thermodynamics with Wilson- 
type quarks, and to study especially the effect of non-zero baryon density. In this report, we will 
highlight the fluctuations at non-zero temperature and density among various topics we are study- 
ing by Wilson quarks. The existence of the endpoint of the first order phase transition in the (r,/^^) 
plane is suggested in phenomenological studies and has attracted much attentions both in theories 
and experiments. Among others, an interesting result has been reported in numerical simulations of 
the quark number susceptibility (the second derivative of the thermodynamic grand canonical po- 
tential Q.) in the framework of the Taylor expansion with respect to jJ-q/T by the Bielefeld-Swansea 
Collaboration [Q]. By calculating the Taylor expansion coefficients of Q. up to 0[{lJ.q/T)^] using the 
p4-improved staggered fermions, they found that the quark number fluctuation increases rapidly as 
jjLq increases in the region near the transition temperature. This suggests indirectly the existence of 
the critical point in the {T,}JLq) plane. It is particularly important to confirm whether the same result 
is obtained using the Wilson fermions which does not resort to the fourth-root trick of the quark 
determinant. Because the odd derivatives of Q. vanishes and the second derivative is the suscepti- 
bility at }JLq = 0, the forth derivative is the leading term necessary to investigate the /x^ -dependence 
of the susceptibility. We calculate the second and fourth derivatives of Q. with respect to pLq, and 
discuss its behavior at finite jXq. 

2. Simulations 

We perform simulations for Nf = 2 at mps/my = 0.80. We adopt the Iwasaki (RG) improved 
gauge action and the clover improved Wilson fermion action: 

^{I5,K,^) = f midctMffe-'^, S, = -I5\ £ colV^v"H^)+ I ^W^f (x) 
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Figure 1: Left: Phase structure of QCD with Wilson-type quarks. Right: Lines of constant mps/my and 
lines of constant T /Tpc- 
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where W^y^{x) and W^^^(x) are 1 x 1 and 1x2 Wilson loops, F^v = (//iv — /^v)/(80' //jv is 
the standard clover-shaped combination of gauge links, j3 = ci = —0.331, co = 1 — 8ci, 

csw = (1 -0.8412j8-')-3/'^, and ^ = jx^a. 

The phase structure of QCD at = with this action has been investigated in Ref. [|l]]. The 
black line (Kc) in Fig. |I| (left) is the chiral limit, on which the pion mass vanishes at zero temper- 
ature. Above this line, a parity-flavor symmetry is spontaneously broken. Numerical simulations 
are performed in the normal phase below Kc. At finite temperature, the parity-flavor broken phase 
becomes narrow, that is the colored region in the upper left of Fig. I (left) for A^, = 4. On the bound- 
ary of this phase, pion mass vanishes. The red line (Kf) is the finite temperature pseudo-critical 
line for Nt = 4, separating the cold phase at small jS and the hot phase at large j3 . 

The relation between the simulation parameters (j8 , K) and the physical parameters is shown 
in Fig. [I] (right). We determine the lines of constant physics (LCP's) by the mass ratio of pseudo- 
scalar meson and vector meson nips/mv, interpolating the data of mps and my at T = in Refs. [|l], 
^ |3|]. The temperature is estimated by the vector meson mass mya using T/my = \/{Ntmva), 
and normalized by Tpc / my, where Tpc is the temperature at the pseudo-critical point on each LCP 
(see Sec. IV in Ref. iQ]). Colored lines in Fig. |l| (right) are the lines of constant mps/mv (LCP's) 
and dashed lines are the lines of constant T /Tpc. {T /Tpc = 0.8 — 2.0.) In this study, we generate 
500 - 600 configurations (5000 - 6000 trajectories) for each simulation point. The lattice size is 
A'^^ X A/, = 16^ X 4. The details of our simulations are given in Ref. [0]. 



3. Taylor expansion of the grand canonical potential 

QCD at finite density is known to have a serious problem called "the sign problem". To avoid 
this problem, we perform a Taylor expansion of physical quantities in terms of }JLq around Hq = 
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and calculate the expansion coefficients, i.e. derivatives of the physical quantities, by numerical 
simulations at }i.q = 0, so that this calculation is free from the sign problem. The calculations 
of the derivatives are basic measurements in the study of QCD thermodynamics, since most of 
thermodynamic quantities are given by the derivatives of partition function ^{T ,}Xu,j^d)-, e.g., 
pressure p and quark number density n are defined by 

where is the chemical potential for the u(d) quark. Quark number susceptibility {Xq) and 
isospin susceptibiUty {Xi) are given by 

Xq _ ( d d \ riu + rid Xi_ _ ( d d \ nu-ng 



Moreover the chiral condensate is defined by the derivative of In ^ with respect to the quark mass. 
We define the Taylor expansion coefficients of the pressure p{T, Hq) for the case Hu = Hd = 

as 

P ^ /^N/Ai«\" i a"hi^ 



(3.3) 



The quark number and isospin susceptibilities for iiu = lid = liq are given by 

^ = 2..+ ,2.,(f)%..., ^jS^^2.-^nc-{tl)\.., ,3.4, 



where 



C„ = 



;U,=0,M/=0 



n\Nl d{^i/Tfd{^q/TY-^ 
The explicit forms of the coefficients are 

C2 = ^=^2, C4 = ^(^4-3^2^), 4 = ^^2, ci = ^(^4-^2=^2), (3.6) 
^2 = m + i&l), i4=(i^4)+4(^3^l)+3(i^|>+6<^2^?> + <^), 

■^2 = (^2) , ^4 = {%) + 2 (^3^1) + <^|> + <^2^?> , (3.7) 
where ^„ =Nf[d"\nAeiM /d^i^] and /i = /i^a. 

^1 =A^ytr[(5M/5jU)M-'], ^2 = A^/[tr[(52M/5Ai^)M-i]-tr[(aM/aAt)M-i(aM/a/i)M-i]],... 
The derivative of the fermion matrix M at = is 

f (1 - y4)f/4(x) - (1 + 74)^/4"^ {x - 4) J for n : odd. 



-K[{\-y,)U,{x) 5^^^^^ + {\ + y,W(.x-A) 5^_^\ forn:even. ^^'^^ 



For the calculation of these operators, the random noise method is used. We generate Nnoise 
of U(l) noise vectors (t]; ^)^^ = T](/,x)5q.^, where T](/,x) is a U(1) random number (t] =e'^;0< 
B < 2k) which satisfies {l/N)Y!i=\ f]{hx)r\*{i,y) = 5x,y for large A'^. a = 1, • • • , 12 is the color and 
spinor index. Then hmA,^<.,(l/Af)£^i L^^ij (^<,a)x,js {^la)yj = ^x,y5p,y, hence 

where is the solution of MX,- « = (• • M"^)?],- «. 



4 



Equation of state for two-flavor QCD with an improved Wilson quark action at non-zero chemical potential 
S. Ejiri 



SB limit 



T/T 



d'(Xq/7^)/d(n/7)' 
d'(x,/r')/d(n/D' 



SB (N =4) 



■ H \ * 



SB limit 



Figure 2: Left: Quark number (black) and isospin (red) susceptibilities. Right: The second derivatives of 
these susceptibilities. 



4. Quark number and isospin susceptibilities 

We calculate the expansion coefficients C2,C4,C2 and c\. The black and red symbols in Fig. ^ 
(left) are the quark number susceptibility Xq/T'^ = 2c2 and isospin susceptibility XijT'^ = 2c\, 
respectively. We also plot preliminary results of the second derivatives of these susceptibilities 
d'^{Xq/T^)/d[fLq/Tf = 24c4 (black symbols) and d^{xi /T^) / d{iXq/Tf = lAc{ (red symbols) 
in Fig. (right). In this study, we choose N„oise = 10 for the calculations of the operators in 



Eq. (3^) except for the operators tr[(5"M/5/x")M^^], where n = \—4. We increase the number of 
noise vectors up to Nnoise = 50 for tr[((9"M/^/i")M^'] to efficiently reduce statistical errors. (See 
discussions below.) 

To check the reliability of the random noise method, we calculate the operators ^„ (?i = 1 — 4) 
using two independent sets of noise vectors with N„oise = 10 on the same configurations . Figure 
|3| (left) shows the time history of the imaginary part of &i and the real part of &2 computed by 
two series of noise sets. The operator is real for n even and purely imaginary for n odd, and 
the expectation value of &\ is zero because an expectation value of imaginary part is always zero 
at /i^ = [^]. It is found that two results of ^2 obtained by different noise sets are consistent 
with each other on each configuration, while two results of ^1 are sensibly different. This means 
that errors from the noise method is dominating in &i with Nnoise = 10. Moreover we found that 
the error from Si\ dominates in C4 and c\ through Eq. ( |3.7| ). Therefore, in reducing the errors for 
second derivatives, it is efficient to increase Nnoise for keeping Nnoise = 10 for other operators. 

As seen in Fig. § (left), Xql^''^ ^^'^ Xi/T''^ increase sharply at Tp^, in accordance with the 
expectation that the fluctuations in the quark-gluon plasma phase are much larger than those in 
the hadron phase. These results agree qualitatively with previous results by staggered-type quarks 
[0, ^. Moreover the result of d'^{xi/T'^)/d{}JLq/TY is quite similar to the results by p4-improved 
staggered fermions [Q]. This suggests that there are no singularities in Xi at non-zero density, as 
discussed in Ref. [||]. On the other hand, we expect a large enhancement in the quark number 
fluctuations near Tpc as approaching the critical endpoint in the {T,}JLq) plane. Although the results 
in Fig. ^ (right) have large statistical errors, the results of d'^{Xq/T'^) / d{}JLq/TY near Tpc show 
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configuration no. 77T 



Figure 3: Left: Time history of 2i\ (top) and SI2 (bottom) obtained by different noise sets at T jTpc = 
0.925, mps/ my = 0.8. Right: Slope of the hne of constant pressure at jj.y — 0. 



quite different behavior from those of Xi- The second derivative of Xq near Tp^. seems to be more 
than three times larger than those at high temperature, suggesting the large fluctuations near the 
critical point. However, the lattice discretization error in the equation of state is known to be large 
for the action in Eq. (11). The short lines in the right end denote the values in the free quark-gluon 
gas (Stefan-Boltzmann) limit, both for Nt =4 and in the continuum. It is clear that we need further 
studies increasing the statistics and decreasing the lattice spacing. 



5. Lines of constant pressure 

It is interesting to compare the line of constant pressure (and energy density) with Tpc{lJLq) and 
the chemical freeze out points [^. Here we consider the pressure constant line in the (T, ^3 ) plane. 



dp 



■^, d{p/T') Ap 
dT T 



AT + 



^dip/ll 
dm 



From this equation, the slope of the /^-constant line at /^^ = is given by 



dr 



d{p/T') 



d{p/T') , Ap 



dT 



+ 



0. (5.1) 



(5.2) 



Using the data oid{p/T'^) /d{iXq/Tf = Xq/T^ in Fig. g p/T'^ and Td{p/T'^)/dT in Ref. [g, we 
estimate the slope of the line of constant pressure at Hq = 0. The results at mps/my = 0.8 are shown 
in Fig. H (right). The slope at /i^ = is about —0.1. This is roughly consistent with the previous 
results by p4-improved staggered fermions at mps/my = 0.7 in Ref. which is denoted by the 
red dot. Moreover the slope of the phase transition line in Ref. is TdTpc/d{iJ.q/T)^ ^ —0.07(3). 
The line of constant p is almost parallel to the phase transition line. On the other hand, the slope 
TAT /&{}jLq/TY of the line of the chemical freeze out in the {T,}JL^) plane are about —0.25 |||]. 
Further studies at small quark mass and large Nt are necessary to compare with experimental results. 
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6. Conclusion 

We reported the cuiTent status of our study of QCD thermodynamics with a Wilson-type quark 
action. The lines of constant physics in the (j8 , K) plane as well as the relation between the pa- 
rameters {^,K) and {T /Tpc,mps/mv) are determined. Simulations are performed on a 16^ x 4 
lattice. The derivatives of pressure with respect to }JLq and /X/ up to 4th order are computed. The 
random noise method is used. For the calculation of 4th order derivatives, the choice of the number 
of noise vector {Nnoise) is important. We discussed the fluctuations of quark number and isospin 
densities. Although the statistical errors are still large, a clear quantitative difference between the 
second derivatives of Xq and Xi is observed. Xq seems to increase rapidly near Tp^ as pLq increases, 
whereas the increase of Xi is not large near Tp^. These behaviors agree with the results obtained by 
p4-improved staggered fermions qualitatively, and with the expectation from the sigma model. 
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